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If the c quark has an anomalous color-electric dipole moment (CEDM), it may serve as a new 
source of CP violation. The strength of such a CP violation depends on the size of the CEDM, 
dc. We propose two effective ways of testing it from the large sample of -0' — >■ J/il) -\- tt"*" -|- tt" at 
the Beijing Spectrometer, and the obtained result, \d'^\ < 3 x 10~^* ecm (95% C.L.), gives the first 
experimentally determined upper bound on the CEDM of the c quark. 
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I. INTRODUCTION 

Searching for new sources of CP violation beyond the 
standard model (SM) is one of the currently interesting 
projects in particle physics. It concerns the explanation 
of the asymmetry between matter and anti-matter in the 
universe. There have been a lot of experimental studies 
on the CP violation in K-meson, B-meson and D-meson 
decays. So far, these experimental results are consistent 
with the SM predictions [1] . 

There have been other possible new CP violation 
sources under consideration, for example, the possible 
electric dipole moments of quarks or leptons [1]. In 
Ref. [2], the CP violation effects in Z boson decays were 
studied. An effective interaction Lagrangian containing 
the relevant CP-violating terms was presented. These 
included the electric and weak dipole moments and the 
color-electric dipole moment (CEDM) of the quarks. In 
the present paper, we arc concerned with the CEDM of 
the c quark. We note that, to the CP-odd correlations 
considered in Ref. [2], this CEDM does not contribute. 
Ref. [3] suggested a test via the decay J/ij] j(j>(f> based 
on a naive quark model calculation. Unfortunately, there 
is no experimental data on the process J/ip — > 70</> so 
far. Ref. [3] only estimated the testing sensitivity from 
the statistics. In this paper, we propose a test via the 
hadronic transition tp' — >■ J/ip + 7r+ + tt^ at the Beijing 
Spectrometer (BES) based on the calculation of QCD 
multipole expansion [4-7] which has proved to be success- 
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ful in many processes [7]. BES has accumulated a lot of 
tp' decays, and the branching ratio for tp' — >• J/'ip+n+TT is 
about 50%, which gives a large sample for testing CEDM 
effect with certain precision providing the first experi- 
mental determination of the CEDM of the c quark. 

The effective interaction Lagrangian including the 
CEDM proposed in Ref. [2] is 

^CEDM = -^<V;ca^'^75 y^cG;;., (1) 

where d'^ is the strength of the CEDM, cr^'' = ^[7'', 7"], 

75 = J7°7^7^7^, Xa is the GcU-Mann matrix for the color 
SU{3), group, and G^^, = d^G^, - d^G^ - gJabcGlGl is 
the field strength of the gluon field. 

^CEDM affects the hadronic transition processes ip' — >■ 
J/tA + 7r+ -t- 7r~ in two folds: 

i. It contributes to the static potential between c and c, 

which causes the mixing between CP-even and CP- 
odd cc bound states, i.e., both ■0' and J contain 
certain CP-odd ingredients such as ipQ'Pi)- 

ii. ^cEDM contributes to the vertices in QCD mul- 

tipole expansion, so that it affects the transition 
amplitudes. 

In this paper, we shall calculate the above two contri- 
butions systematically. 

We first treat Cq^j^j^j as a perturbation to calculate 
its contribution to the c-c static potential with which we 
calculate the energy shifts and CP violating state mix- 
ings. 

We then calculate the theoretical prediction for the 
distribution dr{ip' J/i/j + 7r+ + ■n^)/dMT^T^ and com- 
pare the obtained result with the BES data, which leads 
to an upper bound of d'^. Finally we construct a CP- 
odd operator O from the initial-state and final-state mo- 
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menta in e^e^ — !■ ?/;' — > J/tjj + 7r+ + tt^, and calcu- 
late its expectation value (O) under the amplitude of 
e+e~ ip' J/ip + TT"*" + TT~ . Since the amplitude 
contains a CP-odd piece proportional to d^, (O) is pro- 
portional to d'^. So measuring (O) can provide another 
way of testing d'^. We suggest BESIII to do this measure- 
ment. 

The CEDM interaction Ccedm in (1) is a dim-5 op- 
erator in an effective Lagrangian with a scale parameter 
A ~ TeV beyond which the standard model (SM) should 
be replaced by new physics. The present study is at en- 
ergies far below A and also much below the electroweak 
symmetry breaking scale. See Ref. [8] for a discussion of 
the effective Lagrangian approach in such a case. In this 
paper, we concentrate on studying the contribution of 
^CEDM to the hadronic transition ip' — > J/ipmr. Here 
we would like to explain why other higher dimensional 
CP-odd operators, such as the CP-odd 3-gluonic CP-odd 

operator Og = -^/abcG^.G^^GLe'^''^" [9], need not 
be included in this study. In an effective Lagrangian the- 
ory, an operator with dimension 4 -I- ti is always matched 
by 1/A" from its coefficient. Let us first look at the 
dim-5 operator Ccedm- Comparing it with the dim-4 
SM quark-gluon interaction, we see that the extra di- 
mension of Ccedm comes from the extra derivative on 
the gluon field, i.e., from the gluon momentum k. In 
the transition i// — !• J/tpTrn, k < M^,i — Mj/^ = 590 
MeV. Thus Ccedm is suppressed by fc/A relative to the 
SM quark-gluon interaction. Next we look at the dim-6 
CP-odd operator Og- Comparing it with the dim-4 SM 
triple-gluon interaction, we see that the two extra dimen- 
sions of Og come from two extra derivatives on two gluon 
fields. Thus Og is suppressed by k"^ / A? relative to the 
SM triple-gluon interaction which is of the same order as 
the SM quark-gluon interaction. So, Og is suppressed 
by fc/A < 5.9 X 10""' relative to Ccedm- There have 
been many papers estimating the magnitude of the co- 
efficient C in Og [10] [11], and showing that C is really 
very small. Therefore, theoretically, it is reasonable to 
take only the leading dim-5 CP-odd operator Ccedm 
into account, and ignore all the higher dimensional CP- 
odd operators such as Og in the present study. 

This paper is organized as follows. In Sec. II, we cal- 
culate the contribution of Cq^j^j^.j to the potential be- 
tween heavy quark and anti-quark, and treat it as a per- 
turbation to calculate the energy shifts and state mix- 
ings caused by this contribution. We shall see that both 
J/^p and ip' contain the CP-violating ingr edient ipil^Pi), 
etc. These mixed quarkonium states define the initial- 
and final-state in the transition V'' — > J/i^ + 7r+ + tt~ . 



Then we study the contribution of C^j^jjj^.j, as a new 
vertex, to the QCD multipole expansion amplitudes, and 
calculate all the transition amplitudes up to 0(<ic) in 
Sec. III. In Sec. IV, we calculate the total Mtttt distribu- 
tion (ir(V'' — >■ J/ip + TT^ + TT~)/dM.^Tr and compare it with 
the BES measured result. This leads to an upper bound 
of d'^ which is the strongest bound obtained so far. In 



Sec. V, we propose the alternative way of determining d'^ 
from the experimental data on (O). Sec. VI is a conclud- 
ing remark. 



II. STATIC POTENTIAL AND STATE MIXING 



Derivation of the Potential 



Since d'^ is supposed to be small, the C 



CEDM 



contribu- 
tions to the potential between c and c can be calculated 
by perturbation similar to the derivation of the Coulomb 
potential in quantum electrodynamics [12]. Let the con- 
ventional heavy quark potential be Vq, and the Cq^jjj^j 
contributed potential be Vi. The total potential is 



V = Vo + Vi 



(2) 



In the following, we take Vq to be a QCD motivated po- 
tential, such as the Cornell potential (the simplest one) 
[13] or the Chen-Kuang potential (more QCD, and better 
phenomenological predictions) [14]. Note that the short 
distance behavior of the Cornell potential is the hardest 
(steepest) among the QCD motivated potentials, while 
that of the CK potential is the softest (flattest). Thus 
comparing the results in the two potential models, we 
can see the model dependence of the result. Now we cal- 
culate Vi to lowest order perturbation. The Feynman 
diagrams for the C^^j^j^j contributins to Vi are showm 

in FIG. 1, where the normal vertex is —igs^^-^ for c and 

Wsl'^-^ for c. The shaded circle stands for the CEDVI 
vertex determined by C 



CEDM- 








Plf \ P2 



FIG. 1: Lowest order Feynman diagrams for Vi, where • is 
the normal vetex and the sliaded circle is the CEDM vertex. 



In the static limit, the obtained Vi in the momentum 
representation is 



^1(9) 



A gsd'c {(T -(T)-q 
'3 2 \qP-ie ■ 



(3) 
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Making the Fourier transformation, we finally obtain 



4(cr-o-) •r(5(3)(r) 



4 gs ,/ (<T - <t) ■ r/r 

3 47r r2 



(4) 



The first term serves as a repulsive core, while the second 
term is an attractive force. We shall see later, especially 
from Eq. (AlO) in Appendix A, that to first order pertur- 
bation, the first term does not make any contributions to 
the energy level and the wave function corrections, so that 
only the second term matters. Note that the dimension 
of d'. is m~^. So it is natrual to define 



U — , 



(5) 



nium orbital angular momentum by Al = ±1. This does 
not change the charge conjugation but changes the par- 
ity, i.e., it violates CP. Take the '^S'l quarkonium as an 
example. Vi(r) changes this state to ^Pi. Thus when the 
potential contains Vi{r), the eigen-state is a mixture of 
the ^Si and the ^Pi states. This affects the decays of the 
heavy quarkonia. 



Since Sc is supposed to be small, we can take Vi(r) 
as a perturbation. Let E^i and |n (2''+^)Lj)o (s = 0,1) 
be the energy eigenvalue and the wave function of the 
quarkonium eigenstate with only Vo{r). 



To first order of Sc , the coreection to the energy eigen- 
value is 



where Sc is a dimensionless parameter. Then Vi(r) can 
be expressed by 



3 

4 g^ Sc {a- -a) - r/r 



3 47r TO„ 



(6) 



B. Energy Shift and State Mixing 

We see that Vi (r) contains a factor {a — a) ■ r which 
flips the quarkonium spins by As = ±1 and the quarko- 



Era = Ell + o(n (2^+i)Lj|yi|7i (2^+i)L,7)o. (7) 



We know that |n ^'^^'^'^^ Lj)q is a CP eigenstate, and Vi 
violates CP. So that the diagonal matrix element in Eq. 
(7) vanishes. Therefore there is no energy shift to 0{Sc). 
Energy shift is of 0{S'l). 

The first order wave function correction is 



£0 



xK(2(«±i)+i)(L±l)j)o. (8) 



(±1 



For example, 





|l'5i)o- 


o{l^Pi\Vi\l^Si)o 


-^11 Pi ^l^Si 




12^51)0 + 


o(liPi|l/i|235i)o 

-^23 Si ^11 Pi 


l^Pi) = 


|l'Pi)o + 




^liPi -^l^Si 




|l'i?i)o- 


o{l^Pi\Vi\l^D,)^ 

^ _ 77.0 
^l^Di ^liPi 



3 o(^^S^\V^\l^ , o{l^D^\V,\l'P,), 3 , 

^23Si -^liPi -^l^Di -^liPi 



Here we see explicitly the mixing of the ^Si and ^Pi states. 

The above expressions for the mixed states are not normalized yet. The normalized states are 



(9) 
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ll^Si) = Cl"o |l'5i)o + Cl^ \l'Pi)o + ■■■ 
\2'S^) = C|„" |235i)o + |liPi>o + • • • 

\\^D^) = Cil \\^D^\ + Cl] |liPi)o + • • • , 



(10) 



where 



^10 _ 



,-.20 



o(iiPi|yi|i35i)o 



Si 



^ ^ / o(iiPi|yi|23gi)o 

\ V Si ^ -^11 Pi , 



c. 



11 

20 



i(iiPi|yi|i3gi)o 

^"iPl ^ -^13 Si 



^ / o(l^Pi|^i|l^^i)o 

V -^11 Pi ^ ^I'^Si , 

o(liPi|yi|235i)o 



■^23 Si 



^11 Pi 



,(liPi|yi|235: 



1/0 



p?l 



Pi 



C 



,(iiPi|yi|i35i)o 



Pi 



,(llPl|Fi|235; 



1/0 



-^23 Si -^11 Pi 
I 



o(iiPi|yi|i^iJi)o 

, ^13 Di - K^P^ , 



r'lo _ 



o(l35i|Fi|liPi)o 

-^"iPl ^ ^13 Si 



/ o(liPi|l/i|l^^i)o 

r I ^?3si-^%i , 



^•20 _ 
C-ii — — - 



Ul^Pi|X^i|235i)o 

^ -^23 Si ~ -^11 Pi , 

0(23^11^1 |liPi)o 

^23 Si -^11 Pi 



o(iiPi|yi|i3Di)o 

, ^13^1 - -E^l'lPi , 



\ 



o(liPi|yi|l35i)o 



^i°3si 



P?i 



Pi 



'-'11 — 



o(iiPi|yi|23gi)o \ 

, ^23 Si ^ -^11 Pi / 
0(m|yi|llPl)0 

rO _ po 

^l3P)i -^llPi 



o(iipiiyiim)o 

, -^13131 - ^11 Pi , 



/ o(liPi|l^i|l35i)c 

r I ^i°3si-^?^Pi 



5(liPi|yi|235: 



1/0 



,(llPl|^l|l3Z? 



1/0 



po _ e;0 
-^23 51 ^11 Pi 



^l3p)i -^liPi 



^•12 _ 



\ 



,^11 

'-'12 



o(iiPi|yiim)o 

rO _ po 

^l3_Di -^11 Pi 



,(iipiiFiim)o 

^l3Di -^liPi 



\ 



o(iiPi|yiim)o 

, -^l3_Di - ^11 Pi , 



(11) 



The detailed calculation of the matrix element o(l^^'i|^i|^^^*^^''-^i)o is given in Appendix A. Expand- 
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ing these mixing coefficients up to O^d^/m^) and witfi tfie results given in Eqs. (A15a) and (Al5b), we obtain 
I 



Cil^ = l + 0{St/mi), 



C2"(my,mJ EE C^o^m^, 



-'lO 



10 



3^/3 m, 

8 /2 <5, 



— M "'/ 



^11 



12 



3 V 3 M;^^ - M^p 



c 



12 



m +111 



r 



2 5. 



-"12 



where m^, to^, stands for the magnetic quantum num- 
bers of the initial state orbital angular momentum, the fi- 
nal state orbital angular momentum, and the initial state 
spin, respectively. The numerical values of these mixing 
coefficients can be obtained once a potential model is cho- 
sen. In the CK model [14] and Cornell model [13], the 
values are given in TABLE II in APPENDIX B. 

Actually, even for 5c ~ 0, the states [2^51) and jl^fi) 
are not just the experimentally observed \^') and ['0") 
since the leptonic width of jl'^i'i) is smaller than the 
experimentally measured value by an order of magnitude. 
Usually people believe that the observed [i/;') and \'ip") are 
mixtures of |235i) and [l^Di) [15] [16]: 



W) 



\2^Si)cos6 - |l'''7:>i)sin6' 
12^51) sin 6* 4- [l^L*!) cos 6*. 



(13) 



The mixing angle 9 can be determined by fitting the mea- 
sured leptonic widths. The obtained values of 9 in differ- 
ent models are [15] [16] 



Cornell model : 6 = 10° 

Chen — Kuang model : 9 = 12°. 

Thus the physical and |V'") are [cf. Eq. (10)] 
m^cos9(^Cl° |235i)o + C]i |liPi)o) 



(14) 



+ cos9(cll \l^Di)o + Cl^ \l'Pi)o]- (15) 



3 V 3 m, A/^^-A/^p "/'™.+™. 



So the hadronic transition ip' ^ J / ijj + 7r+ + tt under 
consideration is expressed as 



cos0( Ci^ \2'S,}o + C^^, \l'Pi)o 

sine (cli \l^Di)o + Cll |liPi)o 



E—E^Q 

+ TT + 7r.(16) 



III. CALCULATION OF THE TRANSITION 
RATES AND DETERMINATION OF THE 0{Sl) 
SPA COEFFICIENTS 

A. 0(<5°) and 0{5l) Transitions 

Hadronic transition process is dipicted in FIG. 2, in 
which the transition amplitude contains two factors, 
namely the multipolc gluon emission (MGE) factor and 
the hadronization (H) factor. We shall treat the two fac- 
tors separately. 



h 




FIG. 2: Diagram for a typical hadronic transition. 

We first consider the MGE part. Let and B'' be 
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the color-electric and color-magnetic fields, respectively. 
In the conventional CP-conserving transitions, the MGE 
part contains certain quarkonium-gluon interaction ver- 
tices, e.g., the color-electric dipole (El) interaction, the 
color-magnetic dipole (Ml) interaction, etc. [4-6, 16]: 

El:-d,-E''{X,t), 

Ml : • B''{X,t), 
m^^^-f J{x~X)x ¥ix,t)^j"j^{x,t)d^xll7) 

where and gj^j are effective color-electric and color- 
magnetic coupling strengths, respectively, X is the 
center-of-mass coordinate of the quarkonium, ^(a;,t) is 
the quarkonium wave function at the space-time point 
{x,t). For the CP-violating CEDM interaction vertex. 
Eq. (1) can be written as 

^CEDM = *4V'l7°^°S5y^cGS,-*<^h°a*^75YV^cG, 



The first term is just the interaction between the CEDM 
and the color-electric field. The second term is off di- 
agonal, so that it is nonvanishing only when the lower 

p • er 

components of the quark spinor ±— are taken into 

E + m 

account, i.e., when the quark is moving. So the second 
term means the interaction between the color current and 
the color magnetic field. In the nonrelativistic limit, only 
the first term is novanishing, and the interaction vertex 
is: 



CEDM: -il^{^o- + —a)-E''. (19) 



As in the case of Ml transition given in Ref. [4], after 
certain treatment of the color factor, the effective CEDM 
vertex is 



rrir. 







0\ Aa 

2 



CEDM 



where E = -fE°-. 



2m, 



■{(T-a)-E, 



(20) 



.5c 1^(0 (Ji , --a I na 



(18) 



With these vertices, the transition (16) can be divided 
into the 0((5j?) terms and 0{6l) terms. They are: 



(i) 0{S°) term: 



(a) Ordinary El-El transition of 25* IS with the coefficients cos^CIqC^q: 



E—E21 



+ TT + TT. 



E=Ei, 



(21) 



(b) Ordinary El-El transition of ID IS with the coefficient — sin6'C||C|o- 



— )■ |1 Si)o 



E — i?2i 



E=Ei, 



(22) 



(ii) 0{Sl) terms: 



(c) Ordinary El-Ml transition of 25" IP with the coefficients cos 6'C|^C|^: 



cos 6'C|^ |2^S'i)o 



E—EiQ 



(23) 



(d) Ordinary El-Ml transition of ID IP with the coefficients — sinSCHCig: 



sin6iQ| \VDi 



^ G" |liPi)o 



E — E20 



E—EiQ 



(24) 



(e) Ordinary El — Ml transition of IP IS with the coefficients (cos6'C2o — sin6'C|2)- 
(cos0C]i - sineCl^) |liPi)o \l^Si)o 



E — E2Q 



E^Ei 



(25) 
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(f) Ml-CEDMl transition of 25" IS with the coefficients coa 9C^^^Cl^ 



— !• Cin |1 S'l)o 



E=E2 



E — ElQ 



(26) 



(g) E1-CEDM2 transitions of 25* IS and ID IS with coefficients coseCf^Cl^ and - sin6'Q|Q^, respectively 



cos 6IC|° \2^Si)o 



+ TT + TT. 



E—EiQ 



(27) 



sin0C!2^ |l^i?i)o 



CJq ll"^ 5*1)0 



E — E2 



E=Ei 



(28) 



For a given potential model, there is a systematic way 
of calculating the MGE factors [6] [7]. 

Next we consider the hadronization (H) part, the ma- 
trix elements reflecting the conversion of gluons into light 
hadrons. These matrix elements are at the scale of a few 
hundred MeV, and the calculation is thus highly nonper- 
turbative. So far, there is no reliable way of calculating 
them from the ffi'st principles of QCD, so that we have 
to take certain phenomenological approach. A conven- 
tionally used approach is the soft pion approach (SPA) in 
which the H-factor matrix element is phenomenologically 
expressed in terms of power expansion of the momenta of 
the two pious with unknown coefficients [17]. To lowest 
nonvanishing order, the number of unknown coefficients 
is usually not large, so that they can be determined from 
taking certain experimental input data. This approach 
has proved to be successful in calculating the transition 
rates, the distributions, etc. [6] [7]. However, in the 
present case, it will be difficult if we merely take this 
approach since there are so many kinds of H-factor ma- 
trix elements listed in (21) — (28) containing too many un- 
known coefficients. There are not enough known input 
experimental data to determine them. Another viable 
but cruder approach is the two-giuon approach (2GA) 
proposed and used in Refs. [6] and [7]. In this approach, 
the external pion-fields in the H-factor are approximately 
replaced by two external gluons, so that the matrix ele- 
ments can be easily evaluated as functions of the pion- 
(gluon-) momenta and the two phenomenological cou- 
pling constants and gj^j which can be determined by 
known experimental inputs. Of course this is a crude or- 
der of magnitude estimate. However, it has been shown 
that this crude approach does give right order of mag- 
nitudes of transition rates for many processes [6] [7] [18]. 
The shortcoming of the 2GA is that it cannot give the 
correct angular-dependent distributions such as angular 
distributions, M^^r distributions, etc. This is because 
that the pion is spinless, while the gluon spin is 1. In 
this situation, we shall take both SPA and 2GA in this 
paper for complementarity. 

We can first take the 2GA to calculate all the transi- 
tion rates listed in (21) — (28) which contain two effective 



coupling constants g^ and g^^. Since the total transi- 
tion rate T{ip' J/tjj tttt) is essentially contributed by 
the 0(i5j?) rate which only contains g^, we can take the 
experimental value of T{ip' J/ip tttt) [1], 



r^,, = 304 ± 9 keV, 

B{iP' J/tA tt+tt") = (33.6 ± 0.4)%, 

B(V'' -> J/i^Tr°TT°) = (17.73 ±0.34)%, 



(29) 



as input to determine g^. For the determination of 
gj^j, we can take the 2GA calculated branching ratio 
B{'tp' — > hcTT^) X B{hc rjc'-f) [16] containing g^/g^ to 
compare with the corresponding experimental value [19]. 
For example, in the CK model and the Cornell model, 
the determined g^ and g^ I 9e 



CK model : 



0.523, 



a 



Cornell model : 



5| 
47r 



0.667, 



M 



all 



all 
al 



2.36, 



2.36. (30) 



With the determined g^ and 5^, we can obtain all the 
relative sizes of the transition rates in (21) — (28) with 
only one undetermined parameter (or 8c) left in the 
0{8\) rates. Next, we can take the SPA to calculate all 
the transition rates containing certain unknown coeffi- 
cients. As a perturbation calculation, we ffist obtain 
the 0((5°) transition rate containing three unknown 
coefficients [17] which can be determined by the data 
(29) and M^^ distribution as what is conventionally 
done [20]. Then we calculate the 0((5c) transition 
rates containing several new unknown coefficients. 
Comparing the SPA and 2GA results, we can express 
these unknown coefficients in terms of the known 0((5j?) 
coefficients, so the SPA results of the 0((5^) transition 
amplitudes contain only one undetermined parameter 
d!c (or (5c). We can then calculate the total distribution 
dT{tl.'' — > J/tp + Tr~^ + 7r~)/dAf^7r to compare with the 
BES data for testing d'^ (or Sc)- We shall see that the 
inclusion of the 0{5q) contributions docs improve the 
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fit, and the best fit value of dc is nonvanishing. However, we list the results, 
considering the experimental errors, Sc is still consistent 
with zero. So we can obtain an upper bound on Sc (or 
d'c)- The detailed analysis of this kind of study will be 
given in Sec. IV. 



B. Transition Rates of the 0{S") processes (21) 
and (22) 

The 0(6^) transition rates in Eqs. (21) and (22) have 
been calculated in the pubhshed papers [6] [7] [15]. Here 



The transitions in (21) — (22) belong to the ordinary 
ElEl transitions. The ElEl transition amplitude can be 
written as [6] [7] 



„2 ^ {^p.\r^ \KL){KL\r„^ |$,) 

MeIEI = y. ^ (^^l^-m,-^-mJO) 

I 



where $j {^p) '^^ initial (final) quarkonium state, r 
is the separation vector between the two heavy quarks, f 
is the unit vector of r, r^^ {Tm^) is the component of r in 
the spherical coordinate system with the magnetic quan- 
tum number "T-mi (™m2) ^ {IjO?"!} (cf- Eq. (A3) in 
APPENDIX A), I, s.j.m are, respectively, the orbital an- 
gular momentum, the spin, the total angular momentum. 



and the total magnetic quantum numbers of the quarko- 
nium state, K (L) is the principal (orbital angular mo- 
mentum) quantum number of the intermediate state, and 
EjiL is the energy eigenvalue of the intermediate vibra- 
tional state \KL). The factor {lF^FjF''^F\''^k^i\h^iji''^i) 
can be evaluated using the properties of the spherical 
harmonics [6], and the reduced amplitude f^] ^ ^ is 



f 

•I n 



LP,P^ 



K 



Mi - Ekl 



(32) 



in which i?/, Rp^ and Rkl are radial wave functions of 
the initial, final, and intermediate vibrational states, re- 
spectively. These radial wave functions are calculated 
from the Schrodinger equation with a given potential 

LP P 

model. The values of various /„ in the CK and 

I I F r 

Cornell models are listed in TABLE HI in APPENDIX B. 



In the following, we consider the approaches to the 
hadronization factor (7r7r|i5^£'°|0). 

i. SPA 

In the SPA, the hadronization factor can be generally 
parametrized as [17] 



9l 



llo.{qi)TTp{q2)\E%^,E'L„,^\Q) 



A5, 



■^/2a7|2w^ 



V(2wi)(2w2) 
B 



A 



(33) 



where A,B and C are phenomenological constants, and = {Lu^,q^) is the four momentum of tt^. For a given 
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invariant mass M^^, the A term is angular-independent, the total transition rate are from the A and B terms [20], 
while the B and C terms are angular-dependent [17]. while for |l^S'i)7r7r, the main contribution is 

For |2^5i) |l"^5i)7r7r, the the main contributions to from the C term [6]. Thus the A/^^^ distribution is [6] 



dT{^|J' J/V'tttt) 



SPA 



= \A\^ cos2 



1-41^ 



cos^.^^-(^') 



B . 



B 









fill lj.l\\2 



where 



sm 



-^0 



2 n ^'^J/^ 



2257r3 M. 



c 




A 


1 



J-Q -f (8m^ - A/4^JJ-i H- |/i%(V')r, (34) 



A/, 
2 80 



v/A4\ - Ami, 



i^(Ae + 2ml) + - 4ml) 



l^iM^. + 2Ml^ml + 6mt) + \^{MI^ - 4mJ)(M^. + 6mJ) + (A/^, - 4ml)^ 



in which (A'oi -f^) is the four momentum of the tttt system 



A/^, - A/ 



(35) 



(36) 



The determination of A, B/A and C/^ will be discussed in Sec. IV. 
Integrating over dM^^ in (34), we obtain the transition rate 



SPA 



cos'eG^si^') |/2^o'io(^')l' + sin2 i7c(^') |/i'2"o(^')l' 



(37) 



where 



G 



1 M 



J/V 



~ 87r3 Af , , 



A A/ 



.r2 



1 M 



j/ip 



2m^ 
AM 



\{Ml^ - 2mlfT^ + ^(A/2^ - 2ml)T^ 



B 




A 





dMi 



iml{ml-Kl) + -{Ml,-2ml)^ 



J-o + (8m2-AOJ-^+4^2 , (38) 



in which AAf = M^, - A-f^^^, 

Although this approach can give the reliable M^^ dis- 
tribution, it cannot be applied to other processes in (21)- 
— (28) since other processes contain new unknown pa- 
rameters in their hadronization factors, and there are 
not enough experimental data to determine them. So we 
have to take the help of the 2GA. 

ii. 2GA 

The hadronization factor can be further written as 

N 

where \N) denotes a complete set of intermediate states. 
The 2GA assumes that the color-singlet 2-gluon state 
1-^) = \9'^g^) dominates. In general, the factor (TTTrlg^g^) 



is a function of the pion momenta. Considering the fact 
that the hadronnization factors in most hadronic tran- 
sition processes are all at the few hundred MeV scale, 
the running of (7r7r|(7'^9^) with the momentum in such a 
range of scale is mild. So that we approximately take 
{T^T^\g'^ g'^) ~ const., and the constant can be absorbed 
into the redefinition of the effective coupling constant 
g^. Thus we have 

(7r7r|i?^„^i?^„j0) « (.gV|i?-™,i^-™jO). (40) 

This approximation can be extended to other hadroniza- 
tion factors containing color magnetic field(s). 

The matrix element {g'^g'^\Et^^^E°ira^\^) can be easily 
evaluated: 
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{g^f\Et„,E<L„,^\0) = 



ah 



I 



(41) 



where w (w') and e (e') are the energy and polarization mentary calculations, wc obtain 
vector of the gluon, respectively. After lengthy but ele- 

I 



2 N 2 

9e 



8 (AM)^ 



6 7 277r3 140 



cos^ei/. 



20101 



r 



(42) 



In the following we shall take the 2GA to calculate 
other transition rates in (21) — (28), and compare them 
with (42) to determine their relative sizes. 



C. Transition Rates of the 0{5l) ElMl processes 
(23), (24) and (25) 

These processes contain 5c through the mixing co- 
efficients C20, C^Q, and (cos 06*20 ~ sinflC^j), respec- 



J 



tivcly [cf. Eqs(ll)]. The transition amplitudes are 



123^1) ^ IliPi) 



IliPi) ^ |l35i) 



and ll^Di) \l^Pi) + 7r+ + tt". They ah belong to the 
ordinary ElMl transitions which have been calculated in 
Refs. [6] [16] [18]. Here we present the results as follows. 

For these three processes, let us denote the initial and 
final quarkonium states by and and the spins of c 
and c by and Sg. Then the ElMl transition amplitude 
can be generally written as 



■Meimi = I— 2^ 



12m, 

x{Tr+TT-\E-„,^B-,nJO). 



KLM^ 



M^, - Ekl 



(43) 



Taking the technique of evaluating the spin matrix cl- and rup are the magnetic quantum numbers of the 

ements in the MGE factors (cf. Appendix), and after initial- and final-state, and the definition of fl^'^i'^fi is 
certain lengthy calculations, we obtain i i f f 



(a) \2^Si) |liPi)+7r + 7r 

A/( lo'iQ .1lp^^^ ^ab9E9M/flW , j-OOIn 

6v om, 

X Yj {<ni("^F)am2(mj)}(7r7r|£;_„^B_„jO), (44) 



given in Eq. (32). 

(b) IllPl) ^ \13S^)+7V + 7T 



10-' 



6A/3m, 

{«mi('^F)am2(wj)}(7r7r|£;_„^B_„jO). (45) 



where the definitions of a* and a are given in Eq. (A5), and 
I 



(c) ll^Di) ^ |liPi)+7r + 7r 



Meimi{1'D^ ^ l^Piw) = ^^ii^^(/ii{'^ + f^^l^) ^ {al^{mp)a^,{mj)}{7Tn\E_„^^B^„jO). (46) 

^ C m T m r. 
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Next we take the SPA and 2GA to evaluate the hadronization factor (7r7r|i?_„j^ i3_„j,jO). 
i. SPA 

The hadronization factor {■k'k\E_^_^B_^^\Q) is a second rank pseudo-tensor which is to be expressed in terms of 
the pion momenta in the SPA. To lowest order, the expression is of the form 

9e9m{'^'^\E-„,^B \Q) = /C£;iMie-mi,-m2,-m3 /Trrr^ ' (47) 

where K,eimi ^ phenomenological constant. With this expression for the hadronization factor, the amplitudes 
(44) — (46), after certain evaluation, are 



J^ElMlK^ 6i ^ 1 i^l) - COS6^Cio(j2011 + J20II 



( ^i^. -L ) +<Z2-r)™, (48) 



■^BiMii-L n -J' i iij sni 6^010 (.Jiiio + Jii 



■\/2wi2a7 



plOl \ 

iioJ 



x(-l)'"^-"= ( ) (<Zi-2+g2-i)™, (49) 



-Mi.iMi(l^i^i^l^i^i) = ^g^(cos^^C]o^-sin<i)(/ii°, + /|°,\) 



x(-l)^+"^' ( ) (gi-2 + 92C^i)™, (50) 



where the matrices are the Wigner 3-j symbols. 

Then, after certain calculations, we get the M^^ distributions 





l^Pl)sPA 








5407r3m2 




V'Si)sPA 




















10807r3m2 



where 



|(cos0Cii - sineCl'MZ + /?2°i\)| ^M..-F„J-3, (51) 



4—^ - 3—^ - 1 



f2 ™2 \ _2 / -"-0 



40(M;^^-m^)m^-f -f -1 ■ (52) 



A.r2 I A/2 



Integrating (51) over dM^^, wc obtain the transition rates. For example. 



(53) 



Now the problem is to determine the unknown constant IC^imi ■ there is no accurate enough data to determine 

it, so we should take the help of the 2GA. 

ii. 2GA 

In the 2GA, the hadronization factor can be expressed as 

/IP R ln\^/ IP R ln\ ""'-rn.Wil' X 6^ ( V)) + c^^e^n. ( AQ (? X e(A))-^, 

{TTTrlE^rn^B^rm |0) « {gg\E-rn:,B-rm |0) = /— • (54) 

v2a;v 2w' 
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After certain calculations, we obtain 



iisiA/i(,^ ^ 1 n7r7r)2GA = o ^o_^R 1 777Z~ ) — ?^ — |cos C^C-io (/2011 + /20ii)| S 

2 (AAf)7 



2 (27r)5 \ANmJ 
Compared with Eq. (42), we have the ratio 



AM 



140 



.2GA _ r£iMl(2'5i ^ llPiw) 



BlMl 



r£i£i(235i ^ i-^5i^^) 



f 9m 



= 6 . 

2GA \^''^c9e 



cos^Cjo (/2011 + /2011 



n2 flflll 



2 • 



From (53) and (37) we can get the corresponding ratio 

I^BlA/lP 



T?' 



SPA _ rj;iMl(2'>gl ^ I'PlTTTT) 



ElMl 



r^i^i(235i ^ 13^1 w) 



5407r3m2 



(55) 



(56) 



M /-AM 



■ip' J'lniTT 



SPA 



cos2 G^g(V'') |/iiio(V'')P + sin'e i^yc(V^') Ifl^loi^'W 



(57) 



As has been argued [6, 18], we expect R^mi ~ Reimi- Thus from (56) and (57) we have 

2 



6 



I^BlA/ll 



COS 



2 flflll 
/2OIO 



,2 /) 



1-41^ 



cos^e G^elV/) |/lo"o(^')l' +sin2 l/mo(V'')l' 



1 Af 



540^3 ^^^^^ 



A A/ 



(58) 



So is expressed in terms of „4, S and C appearing in G^g ahd Hq [cf. (38)] which will be determined in Sec. 

I 



D. Transition Rates of the 0{Sl) MlCEDMl processes (26) 

In this transition, 6^ is from the CEDMl vertex. The transition amplitude is 

MmicEDMi = ^^ab E JT^^^ {{^fK^c- S,)rnjKLAQ{KLAQ{s,~ S.UA^j) 

+ {'^f\{Sc - s,),n,KLAQ{KLMj{s^ ~ Sc-)mJ1>/)} (^^|i?_,„,S-„jO) 
= -^^Kbfwio {a™,(m/)<„^(m,) + a^^ (m/)<„^ (m^)} (7r7r|£;_,„,B_,„jO). (59) 

I 

The last step is obtained after certain evaluations of the In (59), the hadronization factor is the same as in the 
MGE factor. ElMl transition. If we take the lowest order SPA expres- 

sion (47), we see that it is antisymmetric in and mj. 
However, the MGE factor in (59) is now symmetric in 

i. SPA 
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and TO2, so that the lowest order SPA does not con- 
tribute to (59). In this case, we should consider the next 
term in the SPA for the hadronization factor in (59), i.e., 



the 0{q'^) term. Considering that {Tnr\E-jn2B-mi\0) is 
a second rank pscudo tensor, we have 



9m x /^^\Tp R in\ _ „■ ^MlCEDMl 
-^(>ab\^^\^-ni2-^-mi\yJ) - « 



CEDMl 



ill - 12) -m^^ 



m-^.—in^. — m 



+ (91 - 92) 'Ilm,q2 



3-i2m^ 



(60) 



r 



Here we have put a factor i reflecting d ^ iq in the This expression is symmetric in and 7712 respecting the 

expansion for convenience, and have introduced a scale Bose symmetry between the two pions. It gives nonvan- 

parameter AM = M^, — Mj^^^ = max(a;i +UJ2) for mak- ishing contribution to (59). 

ing A^j\f dimensionless like other SPA coefficients. With (60), after certain calculations, we obtain 



dV 



MiCEDMii'^^^i ~^ l^- 5'i7r7r)gp^ _ 8 \I^micedmi I ^J/i' 1 /■ooo 12 x- 



K'^{Mt, + Ami) + '^M^M^, - Ami) 



(61) 



\Kl + Ml^-Aml 



-Ml^{Ml^~\mi^)+miKl 



(62) 



Integrating (61) over dM^^, we obtain the transition rate 



''\^MlCEDMl\ ^"c |2 

/ 2010 1 



MlCEDMl[^ ^1 ^I^^JSPA ^ o„3| a i\ rl2 ~~ I/2OI0I / -^4"^"7r7r- 

J2m_ 



37r3|AA/|2 
I 



mt M^, 



(63) 



To determine IC 



the help of the 2GA. 



MlCEDMl 



, we do the calculation with 



ii. 2GA 

The 2GA expression for {TnT\E-m2 B-mi |0) has already 
been given in (54). After certain calculations, we have 



r r'?^'^ ^IS'^ - 32 fgM6c\\.ooo ,2 (AM)^ 

i A/iciJDA/il^ -^1 ^ i 5i7r7rj2GA - '^y'^:^J \J20w\ • 



(64) 



Compared with Eq. (42), we get the ratio 



R'micedmi{'^'^Si -> I'^Sinn) 



^MiCEDMii^''Si^l^Si7rn) 



12 



^EiEii'^^'Si ^ I35i7r7r) 

ffif'^c I./2010P 



2GA 



g%mt 



2 n fill 



cos^ 9f. 



,;„2 n fill 



(65) 
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Taking RI^^cedmi ~ RfncEDMi: determine 

,9m 



MlCEDMl I 



12 



lAMpl/o^p 



5| 



cos^^/iiio ■ -111 



M /-AA/ 



(66) 



in which \I^micedmi\^ expressed in terms of A,B and C appearing in G_a_b a-hd iJc [cf- (38)] 



E. Transition Rates of the 0{Sl) E1CEDM2 processes (27) and (28) 



The two processes (27) and (28) for 2351 l^Siinr and l^Di I^Sitttt belong to the E1CEDM2 transitions (the 
definition of CEDM2 is similar to that of M2 in 77(7r)-transitions [5-7, 18]). The transition amplitude is of the form 

ME1CEDM2 = ^^Sab V —^——{{<i>pTnr\S-E{r-d)\KLM,){KLM,\r'-E\'S>j) 

+ ($^7r7r|r • E\cj,K^j,,){KLMs\S ■ E{r' ■ 9)1$,)} 
= J^^-'^ E ,r ^ F {(j>^|g^,r„,3|XLMM,)(j^£Af,|r:„Jci>,)(^^|(9_,„3i;_„Jii;_„jO) 

+ {'i>p\r^AKLM,){KLM,\S^y,„^\^j){TrTr\E^rn,d-rn,E^rnM ■ (67) 



where S is the total spin operator of c and c. The matrix 
element of S between two quarkonium spin states can be 
evaluated by using Eq. (A20) in the APPENDIX. 

As before, we express the hadronzation factor in SPA 
and 2GA, respectively. 



i. SPA 

In this kind of transition, the hadronization factor is 
a third rank tensor. Considering the Bose symmetry be- 
tween the two pionS; its general form in the SPA is 



5e 
12 



nnlE^mid-msE-rn^] 



9e_ 
12 



TTTtI (9-^3 E-mi )£^-m2 |0) 



■^/2w]2a^ 



92,- 



Here we also put a factor i reflecting d ^ iq for convenience. 
After certain evaluations, we obtain 



d'^ElCEDM2{'^^Si l^SiTnT)gp^ _ 8 |/C^ic'££,^2l^ ^J/i' I fill |2 

U2010 



277r3 



m2 M^, 



+ 2 



(69) 



dr^ic^.m/2(l^^i ^ l-^5iw),P^ _8\JCp,cEDM2\' 61 ^'^.J/i^ \fl^l\'[lM<\F, + F,)-liF, + 22FX (70) 



dM^^ 27007r3 ml M^, 
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where 



1 



J-g - {^2 + 2mlF^ + »7i2 - Kl)Fo} . 



r 



(71) 



Integrating (69) over dM^^, we obtain the transition 
rate of ^(2^51) ^ V(l^S'i)7r7r. 



r_EicBnA/2(2^'5'i — 5- l^S'i7r7r)5p^ 



CEDM2 1 



|2 



J/ip I .111 |2 
^ - - L/2010 



2m_ 



277r3 



(72) 
I 



Next, we determine the unknown constant ^eicedm2 
with the help of the 2GA. 



ii. 2GA 

In the 2GA 



yj2ujT2uj^ 



(7r7r|(a,„^£;^,„j£;^„j0) = 



y^2tJj2w2 



('^l)<^2,mj ('^2)91,-013 + (•^2)^1, -m^ (-^1)92,- 



With this, we obtain 

V r93Q^13c^^^ . ^gg^ I fill 1 2 (AA^)^ 

iiJlCBDA/2(2 Al ^ 1 Al^^)2GA - \-^) I/2OI0I ^Y^. 

Compared with Eq. (42), we have the ratio 



Ti2GA ic^Zq , -i3c^^^ _ ^ E\CEDM2^^^^ ^ 1'"^ 5*17171) 

i^E\CEDM2\^ ^1 ^ i ilTTTTjaGA - ^03 ^ 1 3 Q, 



r^isi(235i ^ 13^1 tttt) 
10 /AA'/y (5g 
27 Ufi j 



2GA 
2010 1 



1^111 ,2 



cos2^/iiiio + 7sin'^/mo 



Taking RWicedm2 ~ Reicedm2: wc have 



I^Bl 



C£;dj\/2I 



5 /AAfV 

4 l^y 



I fill 12 
I./2OIOI 



COS 



:2 fl fill 



20T0 ^ 7 sin^ 0f 
o 



2 ;3 J-111 
1210 



\A\' 



1^1.111 ,2 

7r3 M,„ , 



/2ni 



(73) 



(74) 



(75) 



(76) 



So ^EiCEDM2 expressed in terms of S and C ap- pearing in G^e ^-hd i?c [cf- (38)] 
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Now, in the SPA_^ we have expressed all the unknown 

constants I/C^imiI A^micedm2\ ^^^'^ I^eicedapiI 
occurring in the 0{5l) transitions in terms of A,B 
and C in the 0{6^) transitions. The values of 
^EiMi^^MiCEDM2 and ICej^cedm2 from the values of 
A,B and C obtained from the best fit of the 0{S^) con- 
tribution to the experimental data are listed in TA- 
BLE IV in APPENDIX B. Therefore, the magnitudes of 
the 0{Sl) transitions are characterized by only one pa- 
rameter S^. 

IV. DETERMINATION OF FROM THE BESII 
DATA OF tp' J/ij + 77+ + 77- 

In Sec. Ill, we calculated the transition rates con- 
tributed by the mechanisms (21) — (28) individually, and 
expressed all the unknown coefficients in the SPA in 



terms of the parameters A, i3, and C. Now we are going 
to determine the parameters A, C, and dc from the 
best fit of the theoretical prediction to the experimental 
data. 

The transition amplitudes are functions of the two 
pion momenta qi and q2. For two transition ampli- 
tudes, if their pion momenta dependence belong to the 
same representation of the spacial rotation and reflec- 
tion symmetries, they may have nonvanishing interfer- 
ence term in the M^t-k distribution. So that we should 
take account of such interference terms in calculating the 
Mtt^ distribution. Specifically, there are three kinds of 
interference terms to be considered, namely (i) ElMl- 
CEDM2 interference, (ii) CEDM2-M1CEDM1 interfer- 
ence, and (Hi) interference between the 2^ Si — > I'^S'i tttt 
and l^Di -> l^^i tttt amplitudes in E1CEDM2. So we 
should consider the following total transition amplitude: 



cosdCllC\°MEiEi{'i'Si ^ 1^51) - sxiiOCllCl^^MEiEiiV'Di ^ l^Si) 



C0SeC^^Ci^MMlCEDAIli2'Si ^ I'Si) + COs9CillCillMElCEDAI2i2-'Si ^ l^^i) 

-sm9Cl^Cl°MEicEDM2{l^Di ^ l^Si) + cos9CI^C^^Meimi{2''Si ^ l^Pi) 

coseC^^ - smeCll)cl°MEiMiil'Pi ^ I'^i) - sinOCl^Cl^MEiMiil^Da ^ l^Pi 

I 



.(77) 



We see from Eqs. (34), (51), (61), (69), and (70) that the 
CEDM contributions to dT{ip' J / 4' tttt) / dM^^ are of 
0{5^). So that the mixing coefficients CIq, C|q, and 
Cjl in the first square bracket on the right hand side 



of Eq. (77) should be expanded up to their 0{5^) terms 
which give the CEDM contributions through the normal- 
ization of the mixing coefficients, i.e., we should take 



^10 _ 

I^IQ — 



20 
20 



\ 



\ 



,(iipi|yi|i3gi), 

1 



,(iipi|yi|235i), 
1 



1 |^ll|2 



1 I -,lli2 



(78) 



in the first square bracket on the right hand side of (77). contributions to from in- 



dM^ 
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dividual mechanisms have been given in Eqs. (34), 
(51), (61), (69), and (70). Now we consider the contribu- 



(i) E1M1-E1CEDM2: 



where 



(ii) 



where 



(iii) 



where 



With ah these results, we are ready to determine the 
unknown parameters by the best fit of the theoretical 
prediction to the experimental data. We take the BE- 
ST! data on ^ J/V'tt+tt" [21] based on 1.4 x 10^ 
ij}' events. We proceed the determination by taking the 
CEDM effect as a perturbation. We first take the 0th or- 
der 0{5^) contribution to fit the BESII data. There are 
three unknown parameters A, B and C in the 0th order 
contribution in which A is an overall normalization fac- 
tor irrelevant to the M^^ distribution. It is determined 
by fitting the total transition rate with the experimental 
value [1]. The ratios B/A and C/A do affect the M^^ 



tions from the interference terms. By the same approach 
as in Sec. Ill, wc obtain 



I 

distribution, and they are determined by the best fit of 
the theoretical M^^ distribution with the BESII data. 
The best fit curve in the CK model, as an example, is 
shown by the red dashed line in FIG. 3 together with the 
BESII data. 



The values of |^| determined from (29) and the values 
of B/A and C/„4 determined from the best fit of the 0{5^) 



dV 



E1M1-E1CEDM2 



SPA 



2\/3M, 



457r3m2M^, 



cos 6'/2010-^10 f , 



(79) 



EE V2coseClUnZ + fTii) - V2{coseCll - sin0Cii)(A\% + /mo) " ^^^OClUfH^^ + /f°i\), (80) 

(81) 

(82) 



J-g EE /vo[4-F2 - {2Ml, - ^ml)F, + {Ul^ - ^IK^M, 



J-io ^ ^i^o[16J-2 - {llMl, - 23m2)J-i + (Ia/^ - imlMl, - mlKl)Fo\ 



11CEDM1-E1CEDM2: 



dV 



U1CEDM1-E1CEDM2 



2^/2Af„,„A'/_ 52 



SPA 



J/lp 1- ^ r-^cflcir, flfOOO ^111 77 

Stt^M am '^^'^CEDMll^ElCEDM2 COS » Sm C/2010/1210-^ll ' 

c 



(83) 



-{Kl + Ml,~Aml) 



-MUMl^-AmD+mlKl 



(84) 



i:iCEDM2-ElCEDM2: 

'^^ E1CEDM2-E1CEDM2 



dM^ 



SPA 



8V2M J, <A'U^ A2 
^ J/i> . .2 „ . nflll fill -rr 

~2 \'^E1CEDM2\ COSWSmW/2010/1210-^ 12j 



457r3M^, 



(85) 



EE -{-AMl^F^ + [KlM^ 



TT-K ' TTTT 



AmlK^ - %mlMl,)F, - {Ml, - Ami + ^o)(^ + m^\K\^)F^}. (86) 
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FIG. 3: The best fit theoretical curves of the 0th order contri- 
bution (red dashed fine) and the total contribution (dark solid 
line) in the CK model together with the BESII data on dr(-!/)' — ^ 
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contribution to the M^^ distribution are: 

CK Model : \A\ = 1.641 ± 0.036 

B/A = -0.372 ± 0.006, \C/A\ ^ 0.855 ± 0.020, 
Cornell model : \A\ = 2.09 ± 0.05 

B/A = -0.371 ± 0.006, \C/A\ = 1.000 ± 0.166.(87) 



and K^EiCEDM2 obtained from the best fit values (87) 
are listed in TABLE IV in APPENDIX B. 

As an example, we plot the contributions from vari- 
ous 0((5^) CEDM terms to the M^^ distribution in the 
CK model in FIG. 4 with (5c = 1 (the (5c-independent 
part). We see that the CEDM contribution increases the 
low M^^ distribution and reduces the high M^^ distri- 
bution, which is just the opposite to the 0th order dis- 
tribution. This is the reason why the M^^ distribution 
of ■0' — >■ J/tjjTT^TT~ can sensitively determine the CEDM 
parameter Sc- 

Next we take into account the CEDM contributions 
Eqs.(34), (51), (61), (69), (70),(78),(79), (83), and (85) 
to make the best fit of the total contribution. The best 
fit curve of the total contribution is shown by the dark 
solid line in FIG. 3. The numerical result shows that the 
dark solid curve improves the fit a little bit (with slightly 
reduced value) although the difference between the 
two curves is too small to be visible in FIG. 3. The best 
fit dark solid curve determines the best fit value of Sc 
(d'c) listed in TABLE I in which the error bars of 6c {d'c) 
are determined from the experimental error bars in the 
M^^ distribution. Note that the best fit value of 6c is 
nonvanishing. However, considering the error bars in 
TABLE I, the obtained Sc is still consistent with zero. 
So the model dependence of the present approach is not 
serious. 



The values of the SPA coefficients IC^-^^j^j^. I^micedmi^ 



TABLE I: The Best fit values of Sc (rfc) in the CK model and the Cornell model. 





CK model 




Cornell model 


68% C.L. 


95% C.L. 


68% C.L. 


95% C.L. 


\Sc\ 0.025±0.295 


0.025±0.420 


0.078±0.373 


0.078±0.544 


141 (e-cm) (0.110±1.300)xl0" 


(0.110±1.851)xl0"" 


(0.276±1.320)xl0" 


(0.276±1.926)xl0"" 



For instance, the 95% C.L. upper bound of d'c is 

CK model : |4| < 1.96 x lO"^"* e • cm 

Cornell model : |4| < 2.20 x 10"^"^ e • cm. (88) 

So the model dependence of the present approach is 
roughly 12%. 

We would like to mention that, in Eqs. (58), (66), and 
(76), only the absolute values of the SPA coefficients 
^EiMi^ ^MiCEDMi ^'^'^ /^EiCEDM2 are determined. So 
that there is still an uncertain sign in (79) and (83). 
Actually, If we were able to calculate the hadroniza- 
tion matrix elements from the first principles of QCD, 
there would not be such sign uncertainties. The present 
sign uncertainties are due to the plicnomenological ap- 
proach to the haronization factors taken in this paper as 



lacking of reliable QCD evaluation of the hadronization 
matrix elements. In TABLE I, we only take the sim- 
ple case that all the SPA coefficients are of the same 
sign. Now we consider how will the final result affected 
if they have different signs. First we see from FIG. 4 
that the contributions of (83) is so small that its uncer- 
tain sign only causes negligible effect in the total M^^ 
distribution. Thus only the uncertain sign in (79) mat- 
ters. If we take ICeimiICeicedai2 < in (79), the total 
CEDM contribution to M^^ distribution will be reduced, 
and thus the determined \Sc\ (Mel) ^'^^^ t)e larger. Fortu- 
nately, we see from FIG. 4 that the ElMl-CEDMl con- 
tribution is smaller than the individual CEDM contribu- 
tions. Our calculation shows that the determined \Sc\ at 
68% C.L. wiU change to 0.047 ± 0.383 for the CK model 
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FIG. 4: Various CEDM contributions to the Mtt-k distribution for 
a given Sc in the CK model. 



and 0.123 ±0.511 for the Cornell model. Then the upper 
bound of |c?^| will change to 2.63 x 10~^^ ecm for the CK 
models and 3.09 x 10~^^ ecm, i.e., the uncertainty of the 
upper bound is 34% and 40% for the CK model and the 
Cornell model, respectively. Therefore the effect of the 
uncertain signs is not so serious. We conclude that the 
95% C.L. upper bound of\dl^\ determined from the BESII 
data is 



Kl < 3 X 10-"ecm. 



(89) 



This is the first experimentally determined upper bound 
of the CEDM of the c quark. 

The BES detector has already been updated to BESIII 
with the efficiency of measuring low momentum pions 
significantly improved relative to BESII. So far BESIII 
has accumulated 1.06 x 10^ ip' events, and will be able 
to accumulate 7 x 10* ip' events in the summer of 2012. 
That will be a huge sample. We expect that the new 
BESIII data may determine Sc to a higher precision. 



V. THE CP ODD OPERATOR O 

We can propose another way of determining Sc {d'^) 
linearly from the data of ijj' — > J/ipTnr. Consider the 
process 



e~^e 



(90) 



Let p{—p), qi and q2 be the unit vectors of the momenta 
of the positron (electron), tt"*" and 7r~, respectively. For 
unpolarized e"*" and e~ in the overall cm. system, the 
initial state is then — in the sense of the density matrix 
— CP-even. Therefore any nonzero expectation value of a 



CP-odd correlation of the final state particles is an unam- 
biguous indication of CP violation. With our assumption 
of the CEDM of the c quark, the expectation values of 
the CP-odd operators will be linear in Sc{d'c)- On the 
other hand, to the expectation values of CP-even opera- 
tors, the CEDM can only contribute in 0(6^) or higher 
even powers. We shall now construct a CP-odd operator 
for the reaction (90) following Eqs. (3.20) in Ref. [2] 



0=p-{qi-q2)p- 



qi X q2 

\qi X 92! 



(91) 



Then we define its expectation value which is an experi- 
mental observable: 



^ ^ J PM'Kid^MM'ii'' J/V'71"7r), 



(92) 



where Pj^jij^i is the density matrix, and M (Af) stands 
for the magnetic quantum numbers of "0'- At the "0' res- 
onance, the energy of e"*" (e~) is M^,/2 wihch is much 
larger than the electron mass. Thus the colliding e+ (e~) 
behaves essentially as a massless fermion. With the stan- 
dard couplings for the process e+e" — > 7* — )■ -0', a right- 
handed e"*" can only annihilate with a left-handed e~ and 
vice versa. The resulting density matrix for is 



M',M e {1,0, -1}. (93) 

See Sec. 2.1 in ref. [2] for the analogous process e+e~ — )■ Z 
and set t/^^ = there to obtain Eq. (93). In this case 
the normalization constant N in (92) is just the total 
transition rate r('0' — > J/0 tttt) obtained in Sec. IV. Since 
the CEDM contributions to r(0' — >■ J/0 tttt) is negligibly 
small as can be seen in FIG. 3, we can simply take the 0th 
order transition rate (37) or even the experimental value 
r(V'' ^ J/V-tttt) = 156.04 ±5.78 keV [cf. Eq. (29)] for 
the normalization constant in the following calculation. 

Since O is CP odd, a nonzero (O) can only be con- 
tributed from the CP odd part of dTj^^j^, {ip' J/^p tttt) , 
i.e., the interference terms between the ElEl transition 
amplitude and the CEDM transition amplitudes. From 
the angular part of the phase-space integration, we can 
see that only the l^Di — ^ 1^5*1 tttt part in the ElEl tran- 
sition gives nonvanishing contribution to the ElEl-ElMl 
interference term in (92), while both the 2'^ 5*1 — !• l^Si tttt 
and l^Di — )- l^^i tttt parts in the ElEl transition can 
give nonvanishing contributions to the ElEl-MlCEDMl 
and E1E1-E1CEDM2 interference terms. Thus the result 
will take the form 
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(O) = 



A 



MlCEDMl-^ElEl-MlCEDMl 



1- ^2S^1S-1D^1S 
I^E1CEDM2-^EIEI~EICEDM2 



C_ 



'^ElMl-^ElEl-ElMl I 



-rlD^lS-2S^lS 
MlCEDMl-^ElEl-MlCEDMl 



^El 



-rlD^lS-2S^lS 
CEDM2-^E1E1-EICEDM2 



(94) 



r 



where the I's are the phase-space integrations of the in- 
terference terms which can be calculated from the ap- 
proach similar to those in Sec. IV. So, with the measured 



value of (O), we can determine 5c from (94). Our ob- 
tained results are: 



J 



j2S^lS-2S^lS _ 2 
-^ElEl-MlCEDMl ^ 



r2S-i-lS-lD->lS 



-E1E1-E1CEDM2 



-''ElEl-ElMl 



I 



1D^1S-2S^1S 
ElEl-MlCEDMl 



sin/3(l - cos/3)(gig2 + q^q2) (^q^q^^ + ^^^^^ j duj,duJ2, 
OS yg^^3 ^/20io/i"io J sin/3(l - cosl3){qiql + g^ga) f '7f'72M + -^^1^2 ) dojiduj2, 

- sm —— 

2 9007r3m. 



/1210 / sin/3(l - cos/3)qig2('^i'72 + W2'7i)rfwida;2, 



sin 9 cos 



V J/'i/J Oc „iii 



/i2io/^io / sin/3(3-2cos/3-cos^^)gf(72(9i+'?2)d^idw2, 



'^1D-^1S-2S^1S 
-^E1E1-EICEDM2 



COS P ; 



gi • 92 

|gi||92| 



15757r3 m|? 



I — 



(95) 



where Seimi is given in (80). The complicated Integra- the following numerical results in the the CK and Cor- 
tions can be carried out numerically, and they lead to ncU models: 

I 



CK model : 

A f 

(C) = l^-^^^^^-fl'^-EDA/l - ^■^^^^ElCEDM2 

C 

^A 
Cornell model : 
A f 

(^) = TTvT 1 ^■'^^'^^ MIC ED Ml ~ ^■'^'^'^^E1CEDM2 



O.Q12'i]CEiMl + ^■^'^^•^^MlCEDMl + 0-321/C£;ic£;_DJ\/2 



103 i 
C 

^A 



0.00508/C^iMi + 0.0290/CjvncBDMi + ^■i'^^^EicEDM2 
I 



Sc. 



(96) 



in which the values of |,4| and |C/,4| in the two /C^^j^^^^,^/]^, /C^2C£;r>A/2 ^^'^ §i^^^ TABLEIV in AP- 
models are given in (87), and the values of ^ElMl^ PENDIXB. 
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Now wc come again to the problem of the uncer- 
tain sign similar to those discussed in Sec. IV. We know 
that, in QCD, there is in principle no sign ambiguity 
in the various contributions in (O) in (96). But in 
the present phenomenological approach to the hadroniza- 
tion factors, only the absolute values of the parameters, 

\A\ , \C/A\, |A^EiMlh \^MlCEDMl\j ^^'^ \^E1CEDM2\' 

can be determined as in Sec. III. So, in practice, each 
of the five parameters has an uncertain sign. The un- 
certain sign of A serves as an overall uncertain sign on 
the right-hand-side (R.H.S.) of (96), which makes us un- 
able to determine the sign of Sc- Despite of the overall 
uncertain sign, the values of the curly brackets on the 
R.H.S. of (96) will be affected by the uncertain signs of 
C/A, IC^ij^ji, ^MiCEDMi' ^"^^ ^EiCEDM2- The largest 
term in the curly brackets in (96) is the first term. With- 
out losing generality, we can always take T^micedmi ^ *-* 
with the uncertain sign of A taken into account. If we 

take C / A,K.^ij^^i,K.^i(j^jjj^^2 < 0' curly brackets 
in (96) will take their largest value, 6.43 (CK model) 
and 3.53 (Cornell model). If we take C/A < while 
^EiMi^^EiCEDM2 > 0, the curly brackets will take 
their smallest value, 5.50 (CK Model) and 3.08 (Cor- 
nell model). So the uncertainty of the values of the curly 
brackets caused by the uncertain signs of the SPA co- 
efficients is 25% (CK Model) and 13% (Cornell model). 
This is better than that obtained in Sec. IV. Note that 
the uncertainties in Sec. IV and Sec. V are caused by the 
uncertain signs of different terms. 

Moreover, we may define another related observable. 
Define the asymmetry based on the CP-odd operator O 
as 

. _ N,,,nts{0 > 0) - NeventsjO < 0) 

~ Nevents {O > 0) + N,,,„ts (O < 0) ' ^ ' 

This may also be used to determine dc {d'^) experimen- 
tally. 

So far there is no data on (O). We expect BESIII to 
measure it. 

VI. SUMMARY AND DISCUSSIONS 

If the c-quark has an anomalous color-electric dipole 
moment (CEDM), it will serve as a new source of CP 
violation. In this paper, we study the determination of 
size Sc «) of the CEDM from the BESII data on the 
M^^ distribution in ij/ — > J/ ip tttt within the framework 
of QCD. 

We have first studied the contributions of the CEDM 
to the hadronic transition process ip' — >■ J/ipmr, and de- 
termined the size Sc {d'^) of the CEDM by fitting the the- 
oretical prediction to the BESII experimental data. The 
contributions are in two folds, namely the contribution 
of CEDM to the cc interaction potential which causes 
CP-even and CP-odd states mixing, and the contribu- 
tion of CEDM to the vertices in the hadronic transition 
which affect the Af__ distribution in the transition. Both 



contributions lead to CP violation. Since CP violation is 
supposed to be small, we treat the CEDM effect as a per- 
turbation throughout this paper. We studied these two 
kinds of contributions separately. 

The perturbation calculations of the CEDM contribu- 
tion to the cc potential and state mixings are given in 
Sec. II. The potential is the sum of the conventional po- 
tential Vq and the CEDM contribution For Vq, we 
take two extreme QCD motivated potentials, namely the 
CK potential and the Cornell potential to show the model 
dependence of the present approach. The expression for 
Vi is shown in Eq. (6) in which only the second term 
contributes to state mixings. The obtained normalized 
state- mixing coefficients are given in Eq. (11) [see also 
(12)], and their numerical values are shown in TABLE II 
in APPENDIX B. 

The CEDM contribution to hadronic transition ver- 
tices is more complicated. The 0th order transitions are 
shown in Eqs. (21) and (22) with C^^ = C|5 = = 1. 
The 0{Sl) transitions are shown in Eqs. (23)-(28). The 
transition amplitudes of 0{Sl) transitions are propor- 
tional to Sc, and their transition rates are proportional 
to Sc- To the same order, we must also take account the 
transitions in (21) and (22) with the mixing coefficients 
ofO{Sl) in (78). 

The calculation of the M^^ distribution is quite sub- 
tle. A transition amplitude contains two factors, namely 
the multipole gluon emission (MCE) factor and the 
hadronization (H) factor. For a given potential model, 
there is a systematic way of calculating the MCE factor 
[6] [7] . The calculation of the H factor is a highly nonper- 
turbative problem in QCD. There are two approximation 
methods which can lead to the right order of magnitude 
of the transition rates [6] [7] [15] [18], namely the soft-pion 
approximation (SPA) and the two-gluon approximation 
(2GA). The SPA is a phenomenological approach which 
can correctly describe the angular relation between the 
two pious but it contains unknown constant coefficient (s) 
related to the hadronic matrix element in the H factor. 
The 2GA is a crude approximation which is easy to calcu- 
late but cannot describe the angular relation between the 
two pions correctly. Since we are dealing with the M^^ 
distribution which concerns the angular relation between 
the two pions, we have to take the SPA. However, to 
our experience, the ratios between two transition rates in 
SPA and 2GA are quite close to each other [6] [18]. Thus 
we can use this approximate relation and the 2GA calcu- 
lation to express the 0{Sl) SPA coefficients in terms of 
the 0(5°) SPA Coefficients, cf. Eqs. (58), (66), and (76) 
in Sec. HI. Then we can predict the M^^ distribution by 
treating the CEDM contribution as perturbation. The 
0{Sc) and 0{Sc) contributions to the M^^ distribution 
are given in Eqs. (34), (51), (61), (69), (70). (79), (83), 
and (85). 

We then made a best fit of our 0th order prediction 
(34) to the BESII data (cf. FIG. 3), which determines 
the best fit values of the SPA parameters |^|, B/A and 
\C/A\ shown in eq. (87). Various CEDM contributions to 
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the M^^ distribution arc shown in FIG. 4. We see that 
the behaviors of the CEDM contributions are just the op- 
posite to that of the 0th order contribution (cf. FIG. 3). 
This is why the process ■0' — J irn can sensitively con- 
strain 5c {d'c)- Next we included the CEDM contributions 
to make the best fit up to the 0{S^). It is shown that, 
with the CEDM contribution, the fit is slightly improved 
(with slightly smaller x^), and the best fit values of \5c\ 
and \d'J are liste in TABLE I. We see that the best fit 
value of 1 5c I (l^^cl) is nonvanishing. However, considering 
the experimental errors, it is still consistent with zero. 
The 95% C.L. upper bound of \d'J is shown in Eq. (88) 
which shows that the model dependence of the present 
approach is quite mild. Note that in the present ap- 
proach, only the absolute values of the SPA coefficients 
in the CEDM contribution can be determined. So that 
each SPA coefficient still has an uncertain sign which 
may affect the result. This uncertainty is just due to the 
present phenomenological approach to the hadronization 
matrix element. We have discussed this uncertainty in 
Sec. IV, and the conclusion is that the uncertainty of the 
upper bound is (34 — 40)% which is not so serious. Thus, 
taking this theoretical uncertainty into account, we con- 
clude that the 95% C.L. upper bound of \d'^\ in the present 
approach is jd^l < 3 x 10^^^ ecm [cf. Eq. (89)] which is 
the first experimentally determined upper hound of the 
CEDM of the c quark. 

We have also proposed a second method for determin- 
ing 5c (dc) linearly by introducing a CP-odd operator O 
and measuring its expectation value (O) in Sec.V. We 
have shown in Sec. V that this is a better way of deter- 
mining 5c (dc) experimentally. So far there is no such a 
measurement. We suggest BESIII to do this experiment. 

The state mixings caused by the CEDM of the c quark 
makes the transition rates V'' — ^ he tt" and -0' — > J/V' 7r° 
related to each other. Hence, in principle, the experimen- 
tal data of these two transition rates may give another 
constraint on d^. However, the latest BESIII sample of 
ij}' — >■ /ic7r° based on 106M of -0' events is still rather 
small since the branching ratio of '0' /ic tJ"" is 8 x 10^^ 
[19], i.e., the statistical error in this transition rate is 
significantly larger than that in the present study. Fur- 
thermore, the transition ip' — >■ J/^tt*^ is dominated by 
E1M2 multipole gluon emissions, and the calculation of 
this kind of hadronization matrix element is not so cer- 
tain [22] . Therefore the data of these two transition rates 
cannot provide a strong enough constraint on d'^ compa- 
rable to the one obtained in the present study. So far 
the best experiment for determining the bound on d'^ is 
ip' -> J/V'TTTT at BES. 

The BES detector has already been updated to BESIII 
with the efficiency of measuring low momentum pious 
significantly improved relative to BESII. So far BESIII 
has accumulated 1.06 x 10^ ijj' events, and will be able 
to increase to (7 - 10) x 10*^ V' events in 2012. That will 
be a huge sample. We expect that the new BESIII data 
may determine 5c to a higher precision. 

Estimating the CEDM from some UV theories may be 



interesting for future studies. 
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APPENDIX A: Calculation of the Matrix 
Element o(l^-Pi|Vi|n(2^+^>ii)o 

We show here the explicit expression for the relevant 
matrix element o{l^ Pi\Vi\n^'^'+^\Li)i)o. 

In the nonrelativistic limt, the state |7i(^'*+^)(Li)i)o can 
be decomposed into the radial, angular, and spin factors 



$rUmm. (r, 0,0)= Rnl (r)^," (0, ^'^^1''^'^ , 



(^1) 



where i?„;(r) is the radial wave function obtained from 
solving the Schrodinger equation with a potential without 
Vi, and the spin state x'ff'^^'^ is 



X 



(3) 



(3) 



X1/2X1/2 = 
1 , , 

■y|[Xl/2X-l/2 + X-I/2XI/2J 



1 

71 



(3) 



X 



X-1/2X-1/2 - 

1 r / - - 1 

■y|[x - 1/2X-1/2 - X-1/2X1/2] 



1 

71 



{A2) 



There is a spin-dependent factor (cr ^ cr) ■ r/r in Vi. 
We can take the spherical coordinate 





1 Xl 


+ iX2 
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~ V2 


r 




1 Xl 


- iX2 
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~ V2 
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ro 


^ -fs _ 




r 


r V 






(T± = -^(ctI ± iCT2), CT± = -^(fTl ± 1(72), 

(To = era 5-0 = 0-3 

= ( n n ) ' o-_(ct_) 



cro(o-o) 



\Q 

1 

- 1 




\/2 



(A3) 
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and express it as 



- = (do - CTo) h (ct+ - (J+) — 



r+ 



Finally we evaluate the radial integration. There are 
two terms in Vi{r) [cf. Eq. (6)]. We first look at the first 
term. 



{iao-ao)Yl'ie,^) 

It is easy to see that 
(cr - a)r,iX^^ = 2ar,i{ms)X^ 

am[ms) = {-nis — - — , 1 - 



(A4) 



m 



m — 



-rris ^ ), 

m — — 



H. . ^ , l + m., 
amims) = {-ms ,1 



\ms\,~ms ), 



(A5) 



m = + m = m 
and am{ms), 0^(771^) satisfies 



(AG) 



mjUs 

Thus 



' — {2ao{m,)Yl" 



+2a+im,)Y-\e,<j))-2a-ims)Y^\e,(j))}X^ 

{a -a)- -X' = J— {2aS(™.)r°(0, 0) 
r V 3 

+2a;(7n,)rf 1(0,0) - 2a*_(m,)yii(0,0)} Xi^.{A7) 

This explicitly shows that (cr — a) ■ r/r flips the quarko- 
nium spin. 

Next we evaluate the angular integration. What we 
need to evaluate is the integration of the product of 3 
spherical harmonics. According to the property of the 
spherical harmonics, we have 



Y, 



]';'ie,<t>)Yr{o,m;"''io,^) dn 



If 1 k\ ( If 1 k 



/ \"^/ ™ — TOj 
The values of some relevant 3 — j symbils are 



(A8) 



1 1 

ruf m {) 

1 1 2 





-1 
2 

l"5' 



a/3 



/ (2 + 771/ + m)!(2 - m/ - to)! 
^ V 30(1 + TOf)!(l - TOf)!(l + to)!(1 - to)! 



.(A9) 



8g.(r) 
3 An m. 



0. (AlO) 



Here we have considered the running of the QCD cou- 
pling constant g^ir) in the radial integration. Note that 
gg{r) is governed by asymptotic freedom as r — 0, i.e. 

g'lir) ^^^^ (cf. Eqs.(A14) and (A13) below). 

Eq. (AlO) shows that the first term in Eq. (6) actually 
does not contribute. We should only take into account 
the contribution of the second term in Eq. (6) to the ma- 
trix element. The radial integration from the second term 
contribution is 



3 TO, 



9 sir) 



in 



Rnir)Rn,uir) dr 



3 TO, 



Rl^ir)R^,(r) dr. 



(All) 



The radial wave function R{r) is to be obtained by solv- 
ing the Schrodinger equation. 

To include nonperturbative contributions to g^ (r) near 
the J/'ip and -0' scales phenomenologically, we take the 
CK potential model [14] which has both a clear QCD 
interpretation and successful phenomenological predic- 
tions. The CK potential reads [14] 



y(r) 



16n 1 
+kr, 



2jE 



53 
75 



fir) 



462 ln/(r) 
625 /(r) 



{A12) 



where k = 0.1491 GeV is the string tension related to 
the Regge slope [14], je is the Euler constant, and /(r) 
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/(r) = In 



1 



A- 



MS ^ 



-4.62- (1 



4 A:?— 
MS- 



1 



-cxp<^ - 



15 3 



A 



MS 



A- 



MS 



1 2 



-1 IAms 



(A13) 



in which A-^ = 180 MeV. The nonperturbativc effects 
resides in the phenomenological function /(r). Writ- 
ing the potential (A12) in the standard form V{r) = 
4 a (r) 

+ kr, we read from (A12) that 

3 r 



.9s (^) 
47r 

127r 1 



2jE 



53 
75 



fir 



462 ln/(r) 
625 /(r) 



(A14) 



This running formula will be used in the calculation of 
in (All). 

Putting all the above results together, we obtain the 
expressions for o(l^Pi l^i |n?S'i)o and 0{l^Pi\Vi\l^Di)o: 

o{l^Pi\Vi\nfSi)o imf,ms) = -- ^ 



xl, 



" S 



(A15a) 



3V 3m, 



(A156) 



With these two matrix elements calculated, we can obtain 
all the mixing coefficients in (11) [cf. Eq. (12)]. 

For the spin-1 states. The conventional Cartesian co- 
ordinate representation of the spinal operators are: 



/ 

5*1 = - i 

\0 i 

i.e., 



'0 

52 = I I , 5., = 

-i 0, 



(Si. 



The eigenvectors of, for example, 5*3 are: 




(AU) 




rrif, = — 1 

Xi,„. = i(l-H Xi,„. = ^|0|, (A18) 



2 








where the column is ordered according to i = 1, 2, 3 from top to bottom. 

In the polar coordinate system, we should make the linear combination for the component index (xim )±i 
[(Xim )i i *(Xim )2]/>/2- Thus wc obtain the following eigenvectors in the polar coordinate system 

mg ~ 





where the column is ordered according to to = -1-1, 0, — 1 
from top to bottom. Compared with Eq. (A5), we see 
that 

1 + i 

(XimJ™ = —^a^ims) = Na„rims), (A19) 

where N = {1 + i)/V2 is the normalization factor, and 
N*N = 1. 



Now the term {x*lrnJJ{S^)jk{XlmJk^^ ^he 
E1CEDM2 transition amplitude can be evaluated as 
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^ — mi, — m2,— ms (Xlr; 



ima2 ) 



(A19) 

= e 



mi , — m2 m3 ^m2 (^s2)^m3 (^s3 )-^mi 

I (a(ms2) X a(ms3)) • (A20) 



APPENDIX B: Numerical Results of f^!^,'^'', , 

^ !*■ F'' F 

C:;//, , and x;'s 

Here we list the numerical values of the the mxing co- 
efficients C^ji^", reduced amplitudes f^^in^i : ^'^id the 
SPA hadronization factor coefficients /C's in the CK and 
Cornell models. 

TABLE II: Values of the mixing coefficients C/'" in the CK 
and Cornell models. 



CK model 



Cornell model 



""TTrr 
^11 

^11 
'^12 



0.1165(5c 
-0.07119(5c 

0.1597(5c 
-0.01863(5c 
-0.03906<5e 



0.1U2Sc 
-0.072095c 

0.1669<5c 
-0.018245c 
-0.04091<5c 



TABLE III: Values of the reduced amplitudes i' i 
the CK and Cornell models. 



CK model 



Cornell model 



T 

) 

(GeV-3) 



/2cm) (GeV 
f^Z (GeV-^) 

rill 

Jl210 

(GeV-2) 

f?Z (GeV-2) 
/•no 
/20l: 
/•OOl 
/2011 

/i^2\°i (GeV-^) 
f?^h (GeV-2) 



1 (GeV- 
(GeV- 



8.8715 
0.3869 
-11.2507 
-3.31446 
-6.8977 
4.8399 
4.5643 
-5.8859 
-4.1087 



6.9330 
0.3709 
-8.7701 
-2.9040 
-6.1708 
4.2801 
4.0379 
-5.1934 
-3.6324 



TABLE IV: Values of the SPA hadronization factor coef- 
ficients \K,EiMi\. [ICmicedmiI, and \ICeicedm2\ obtained 
from the best fit values of |^|, B/A, and C/A in the CK and 
Cornell models. 





CK model 


Cornell model 


I^BIMI 1 


3.289 


4.208 


I^A/lCEDA/1 


2.318 


2.625 


\^E1CEDM2\ 


0.532 


0.602 
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